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Abstract It is only now, with low-frequency radio telescopes, long exposures 
with high-resolution X-ray satellites and 7-ray telescopes, that we are be- 
ginning to learn about the physics in the periphery of galaxy clusters. In 
the coming years, Sunyaev-Zeldovich telescopes are going to deliver further 
great insights into the plasma physics of these special regions in the Universe. 

["tI . The last years have already shown tremendous progress with detections of 

shocks, estimates of magnetic field strengths and constraints on the parti- 
cle acceleration efficiency. X-ray observations have revealed shock fronts in 

"^ ' cluster outskirts which have allowed inferences about the microphysical struc- 

ture of shocks fronts in such extreme environments. The best indications for 
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Q ■ magnetic fields and relativistic particles in cluster outskirts come from obser- 

H ' vations of so-called radio relics, which are megaparsec-sized regions of radio 

C/3 , emission from the edges of galaxy clusters. As these are difficult to detect 

^1 ' due to their low surface brightness, only few of these objects are known. But 

they have provided unprecedented evidence for the acceleration of relativistic 
particles at shock fronts and the existence of /zG strength fields as far out as 
^ ' the virial radius of clusters. In this review we summarise the observational 
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1 Introduction 

Clusters of galaxies are filled with a dilute (rig < 0.1 cm"^), hot (T > lO^K) 
gas, the so-called intracluster medium (ICM). In this review, we are primarily 
concerned with regions outside of cluster cores, close to the virial radius (^1 
Mpc). The ICM contains the majority of baryons in galaxy clusters and has 
a temperature close to the virial temperature of the gravitational potential 
of the cluster. 

Shocks play a fundamental role in the evolution of the ICM and affect 
clusters in two ways: First, shocks thermalize the incoming gas allowing it 
to virialise and providing much of the pressure support in baryons. Secondly, 
because these shocks are coUisionless features whose interactions in the hot 
plasma are mediated by electromagnetic fields, it is possible for a fraction of 
the thermal distribution of particles to be accelerated and transformed into 
nonthermal populations of cosmic rays (CRs) through diffusive shock accel- 
eration. This process results in a part of the kinetic energy to be converted to 
non-thermal components. Non-thermal populations of particles can in prin- 
ciple be observed in the hard X-rays through inverse Compton scattering 
of the Cosmic Microwave background, in the gamma-rays through hadronic 
interactions of non-thermal protons with the thermal component and in the 
radio band through synchrotron radiation of the electrons. The latter has 
been observed in galaxy clusters in the form of diffuse radio sources. These 
sources are usually su bdivided into two classes, denoted as 'radio haloes' and 
'radio relics' (see e.g. iFerrari et al.l (|2(308l for a review). 

Radio haloes are diffuse, low surface brightness (~ 10~^ Jy/arcsec^ at 
1.4 GHz), steep-spectrurro (« > 1) sources, permeating the central regions 
of clusters, produced by synchrotron radiation of relativistic electrons with 
energies of ~ 10 GeV in magnetic fields with B ~ 0.5 — 1 ^G. Radio haloes 
represent the best evidence for the presence of large-scale magnetic fields 
and relativistic particles throughout the intra-cluster medium (ICM). They 
are mostly unpolarized and have diffuse morphologies that are similar to 
those of the thermal X-ray emission of the cluster gas. Their origin and 
the r ole that shocks play in their formation are still subject to debate, (see 
e.g., iPetrosian and Bvkov|[2008i iPetrosian et al.|[2008l iBrunetti et al.ll2009l 
for recent review). Unlike haloes, radio relics are typically located near the 
periphery of the cluster. They often exhibit sharp emission edges and many 
of them show strong radio polarization. For some relics there is now good 
evidence that these electrons are produced by shock waves and we will present 
this evidence in Sec. 4. 

If a significant fraction of the pressure of the ICM is provided by CRs, this 
will affect gas mass fraction estimates. Since gas mass fractions in clusters 



are used to make inferences about dark energy, a knowledge about the non- 
thermal content of galaxy clusters is important for their use in precision 
cosmology (e.g. Vikhlinin et al. 2010). 

There are three types of phenomena that create shocks in the ICM. The 
first class of shocks are accretion shocks: At very large off-center distances, 
typically several Mpc, cosmological simulations predict that intergalactic 
medium (IGM) continues to accrete onto the clusters through a system of 
shocks that separate the IGM from the hot, mostly virialized inner regions. 
The IGM is much cooler than the ICM, so these accretio n (or infall) shocks 
should be strong, with Mach iiumbe rs M ~ 10—100 fe.g.. lBertschingeiill984 
IZeldovichlll972l iMiniati et al.l 120001 ). As such, they were suggested as sites 
of effective cosmic ray acceleration, with consequences for the cluster energy 
budget and the cosmic gamma-ray background. 

The second class of shocks are merger shocks: as subhalos fall into the 
main clusters, they create moderate-strength shocks that allow the gas to 
virialise. If an infalling subcluster has a deep enough gravitational potential 
to retain at least some of its gas when it enters the dense. X-ray bright region 
of the cluster into which it is falling, a merger shock is produced. In Sec. 2, 
we will discuss some prominent observations of this type of shocks. 

The third class of shocks are AGN shocks: Powerful AGN are observed to 
inflate large bubbles in the IC M, which may generate shocks within the cen - 
tral few hundred kiloparsecs (JBrtJggen et al.l I2007L ISimionescu et al.l l2009t) . 



Like merger shocks, these shocks have moderate Mach numbers, and have 
been detected in cluster centres, such as in Virgo and Hydra A. These shocks 
are very important because they convert the energy produced by the central 
supermassive black hole into thermal energy and, thus, con iplete the feed- 
back loop which is regarded as essential for galaxy evolution ()Cattaneo et al.l 
|2009[ ). As these shocks are primarily found in cluster centres, they will not 
be the subject of this review. 

The setup of this review is as follows: In Sec. 2, we introduce some phe- 
nomena that we will cover in this review and explain some basic terms. In 
Sec. 3, we describe our current knowledge of cluster shocks from X-ray and 
7-ray observations. Then we summarise what we can learn from radio obser- 
vations in Sec. 4, including estimates of magnetic fields in cluster outskirts 
from recent observations of radio relics. In Sec. 5, we review recent results 
concerning structure formation shocks from cosmological simulations. In Sec. 
6, we then discuss the origin of magnetic fields and relativistic particles in 
shocks. Finally, in Sec. 7, we summarise insights gained in the last years on 
the physics of shocks in such dilute environments as found in cluster outskirts 
and identify open questions. 



2 Fundamentals 

Shock structure: Shock waves are usually regarded as a sharp transition 
between a supersonic (and super-Alfvenic) upstream fiow and a subsonic 
downstream fiow. In a rarefied hot cosmic plasma the Coulomb collisions are 
not sufficient to provide the viscous dissipation of the incoming fiow, and col- 
lective effects due to the plasma flow instabilities play a major role. Hence, 



cosmological shocks in a rarefied liiglrly ionised plasma are collisionless. Pro- 
tons, whose thermal velocity is lower than the shock velocity, are heated 
dissipatively at the shock layer. However, the equilibration of electrons may 
differ substantially from that of the ions because of their smaller gyroradii. 
The amount of collisionless heating of electrons depends o n the shock veloc- 
ity an d the Mach number (see e.g. ISchwartz et al. 1988 , Bvkov and Uvarovl 
Il999[ ). Then the Coulomb equilibration eventually occurs on scales that are 
ord ers of magnitude longer than the collisi onless heating scale length (see 
e.g. lFox and Loeblll997L iBvkov et al.ll2008b[) . Under such circumstances, the 
shock structure can be more complex because electrons and ions can remain 
out of equilibrium for a while. Then such a shock consists of an ion vis- 
cous jump and an electron-ion thermal relaxation zone. We will revisit the 
microphysics of cluster shocks in Sec. 6. 

Shock acceleration: Shocks can also change the momentum distribution of 
a gas by accelerating a fraction of the particles to very high energies. This has 
been observed in astrophysical shocks on a wide range of scales. Examples 
are the shock waves that form when the solar wind collides with planetary 
magnetospheres, shocks surrounding supernova remnants in the interstellar 
medium and shocks in active, extragalactic sources such as quasars and radio 
galaxies. The resulting energies reach up to about 100 keV in the shock in the 
Earth's magnetosphere to about 100 TeV in supernova remnants. Particles 
gain energy whenever they cross the shock front in either direction. Scattering 
processes make the velocity distribution isotropic with respect to the local 
rest frame on either side of the shock front. Consequently, there is a certain 
probability that a particle that has just gone through the shock and been 
accelerated will diffuse back through the shock to where it came from. Every 
time the particle crosses the shock it will gain energy. Moreover, the increase 
in energy is the same no matter from which direction the particle crosses the 
(non-relativistic) shock. Particle passing back and forth through the shock 
can thus attain very high energies. This process is called first-order Fermi 



(Fermi-I) accelerat i on or diffus iv e shock acceleration (PSA: iKrvmskii 
Axford et al."l977', 'Bell"l9 78allbl. iBla ndford and Ostrikei^'l978 l iDrur" 



1977 . 



1983, 



Blandford and Eichler 1987r iJones an d Elhson 1991, Malkov a nd O'C Drurvl 
200 1[ ). The cosmic rays thus generated can penetrate far into the upstream 



flow, thus creating a shock precursor. The particles decelerate the flow and 
preheat the gas. If the cosmic ray pressure inside the shock becomes compa- 
rable to the thermal pressure, the cosmic rays will also modify the structure 
of the shock. Finally, shocks can efficiently generate strong fluctuating r nag- 
netic fields in the upstre am region due to cosmic-ray driven instabilities OBelll 
l2004llBvkovet al.ll201lh . 

Cluster magnetic fiel ds: The origin of m agnetic fields in galaxy clusters 
remains unclear (see e.g. iDolag et al.l ( 20081 ) for a review). It has been sug- 



gested that they are primordial; i.e., a seed field that has formed prior to 
recombination is subsequently amplified by compression and turbulence. Al- 
ternati vely, it has been proposed that the magnetic field is of proto-galactic 



origin (jKulsrud et al. 1997|, iRyu et al.. i2008) or that it had been produced 



by a cluster dynamo (JRoettieer et al.l Il999t ). Others conclude that mag- 



netic fields can be pro duced efficiently in shocks by the Weibel instability 
(iMedvedev et al.l2004[) or by small-scale plasma instabilities ( Schekochihin et ahl 



|2005[ ). Finally, it has been suggested that the intergalactic medium has been 
magnetis ed by bubbles of r adio plasma that are ejected by active galactic 
nuclei (iFurlanetto and Loeb 2001 1. Cluster magnetic fields are the subject of 
a separate chapter in this review and in this review we will focus on magnetic 
fields in cluster outskirts only. 

Turbulence: Mergers, the motions of galaxies and the action of AGN pro- 
duce turbulence in the ICM. Turbulence may have a deep impact o n the 
physics of the ICM (see e.g. ISchekochihin et ahl l2005l iLazarianl l200(J l and 
on the properties of non-th ermal components in galaxy clusters (see e.g. 
iBrunetti a nd Laza rian|[201l[ ). The presence of turbulent gas motions in the 
ICM is suggested by measurements of the Faraday Rotation of the polar- 
ization angle of the synchrotron emission from cluster radio galaxies. These 
studies show that the ma gnetic field in the ICM is tangled on a broad range 
of spatial scales (see e.g. IVogt and Enfjlinl l2005l ) suggesting the presence of 
super-Alfvenic motions in the medium. Also independent attempts from X- 
ray observations of a number of nearby clusters, based on pseudo-pressure 
maps of cluster cores and on the lack of evidence for resonant scattering ef- 
fects in the X-ray sp ectra, provide hints of turbulence in the ICM (see e.g. 
IChurazov et al.ir2004[) . Important constraints on the fraction of the turbulent 
and thermal energy in the cores of clusters are based on the analysis of line 
broadening in the emitted X-ray spectra of cool core clusters ( Sanders et al.l 



[2010). In the outskirts of clusters, simulations predict that the tu rbulent pres- 



sure s upport increases relative to the therm al pressure support ( Vazza et al.l 
l2ninllPaul et al.ll2ni(llNagai and Ladl2nilh . 



3 Cluster shocks in X- and 7-ray emission 

3.1 Cluster shocks in X-rays 

Regions of high-entropy gas in clusters have been observed by ROSAT, 
ASCA, XMM-Newton and Chandra and have been interpreted as the re- 
sult of shock heating (e.g., Bclsolc et al.ll2005[) . However, only a handful of 
shock fronts, exhibiting both a sharp gas density edge and an unambigu- 
ous temperature jump, have been found : in the " bullet cluster" , IE 065756 
(JMarkevitch et al. 2002), Abell 520 (Ma rkevitch et al. 2005]), and AbeU 2146 
(^Russell et al.ll2010l ). Such discoveries are so rare because one has to observe 



a merger when the shock has not yet moved to the low-brightness outskirts, 
and is propagating nearly in the plane of the sky, to give us a clear view of the 
shock discontinuity. In addition, the shocks in A520 and IE 065756 have Mach 
numbers of M ^ 2 — 3, which provide a big enough gas density jump to enable 
accurate deproj ection. Merger shock fro nts have been found in other c lusters, 
e.g., in A3667. (iFinoguenov et al.ll20"Toh . A754, (JKrivonos et al.ll2003^ ■ A521, 
(JGiacintucci et al.ll2008l) " 



Recently, iFinoguenov et aP ( 2010t ) analysed XMM-Newton and Suzaku 



pointings of the cluster Abell 3667, including the NW relic (see Fig. 1). They 
find a sharp X-ray surface brightness discontinuity at the outer edge of the 
radio relic, and a significant drop in the hardness of the X-ray emission at the 
same location. This discontinuity is consistent with a Mach number 2 shock 
moving at 1200 km/s associated with the merger. Kinetic energy is being 
dissipated in this shock at a rate of dEKE/dt ^ 1.8 x 10^^ erg/s. The total 
non-thermal luminosity of the relic is Lnt ~ 3.8 x 10'*^[(3.6 ^G/B)'^ + 1] erg 
s^^. If this energy is provided by the shock acceleration of electrons, then a 
fraction e ~ 0.2 % of the kinetic energy fiux through the shock is converted 
to the non-thermal luminosity. Shock acceleration at the outer edge of the 
relic and radiative losses as the electrons are advected away from the shock 
can explain the rapid steepening of the radio spectral index with distance 
from this edge of the relic. Alternatively, the surface brightness discontinuity 
and hardening of the X-ray spectrum might be due to Inverse Compton (IC) 
emission from the relic. In this case, the magnetic field in the relic is about 
3 fj,G. Since the observed X-ray excess from the relic is an upper limit to 
IC emission, this yields a lower limit on the relic magnetic field of > 3 //G. 
This is a remarkably strong magnetic field at this large projected distance 
(2.2 Mpc) from the cluster center, but is consistent with Faraday rotation 
through the relic observed towards two background radio galaxies. 

The bow shock in IE 065756 offers a unique experimental setup to deter- 
mine how long it takes for post-shock electrons to come into thermal equi- 
librium with protons in the ICM (Markevitch et al. 2006). Currently, X-ray 
observations only yield the electron temperature but not the ion temper- 
ature. However, one can use the measured gas density jump at the shock 
front and the pre-shock electron temperature to predict the post-shock adi- 
abatic and instant-equilibration electron temperatures, using the adiabatic 
and the Rankine-Hugoniot jump conditions, respectively, and compare them 
with the data. Furthermore, if the downstream velocity of the shocked gas 
fiowing away from the shock is known, one can deduce how the fiow spreads 
out the electron temperature along the spatial coordinate in the plane of the 
sky. The Mach number of the IE 065756 shock is sufficiently high, such that 
the adiabatic and shock electron temperatures are sufficiently different (for 
M < 2, they become close and difficult to distinguish, given the temperature 
uncertainties). Moreover, the distance traveled by the post-shock gas during 
the time, is well-resolved by Chandra. Markevitch et al. (2006) found that 
the temperatures are consistent with instant heating; equilibration on the 
coUisional timescale is excluded, although with a relatively low 95% confi- 
dence. The equilibration timescale should be at least 5 times shorter than 
the Coulomb time. We will revisit the theory of dissipation in coUisionless 
shocks in Sec. 6. 

The low surface brightn ess of the cluster outskirts is challenging for cur- 
rent X-ray telescopes (e.g. iBohringer and Werneij[2010l lEttori and Moleiidil 
|2010[ ). Projected temperature profiles to ~ 0.5r2oo of 1 5 nearby (z < 0.2) 
clusters were derived from XMM-Newton observations bv lPratt et al.l ( 2007t ) 
(see also I Vikhlinin et al.ii2005i) . 




Fig. 1 Wavelet reconstruction of the 0.5-2 keV X-ray mosaic of A3667 from XMM- 
Newton. There is a sharp discontinuity in the X-ray surface brightness to the NW; 
the blue ellipses show the location, name and size of three main galaxy components 
(KMM5, KMM2, and KMM4) of A3667 from the analysis of Owers et al. (2009). 
The black cross marks the center of KMM2 group. The white contours show the 
wavelet reconstruction of the SUMMS 843 MHz radio image in the region of radio 
relics. The contours are drawn at the 1, 3, 9, and 20 mjy/ beam levels. From 
Finoguenov et al. (2010) with kind permission. 



The low and well-constrained background levels make Suzaku the prime 
instrument for the study of X-ray emission from cluster outskirts. Suzaku is 
in a low orbit within Earth's magnetopause providing significantly lower and 
stable particle background compared to Chandra and XMM-Newton. How- 
ever, the calibration of the instruments can sometimes be diffic ult. The low 
surfac e brightness outskirts of A2204 around r2oo were studied bv lReiprich et al.] 
(|2009[ ) with Suzaku. iSimionescu et a l. (2 ID determined ICM properties out 



to the edge of the Perseus Cluster. Contrary to earher findmgs, the cluster 
baryon fraction is consistent with the expected universal value at half of the 
virial radius. The apparent baryon fraction exceeds the cosmic mean at larger 
radii, suggesting a clumpy distribut ion of the gas for rad ii larger than half of 
the virial radius. Also with Suzaku, George et al. (2009) studied X-ray emis- 
sion from the outskirts of the cluster of galaxies PKS 0745-191 at z = 0.1028. 
They determined radial profiles of density, temperature, entropy, gas fraction, 
and mass beyond the virial radius out to ^ 1.5r2oo- The temperature is found 
to decrease by roughly 70 % from 0.3-1 r2oo- Suzaku observations provide evi- 
dence for departure from hydrostatic equilibrium around and ev en before the 
virial radius as it is seen in the galaxy cluster Abell 1795 (e.g. iBautz et al.l 
|2009[ ). Suzaku observations of PKS 0745-191 revealed one shock front possibly 
indicating that this cluster is still accreting material, likely along a filament. 
Evidence for nonthermal pressure support suggests that bulk motions from 
merger activity could be mak ing a significant con tribution to the gas energy 
in the outskirts of this cluster ([George et al.|[2009[) . There is some evidence for 
the presence of an excess above the dominating thermal e mission of the ICM 
hot gas in both X-rays and EUV (JRephaeh et al.ll2008l iDurret et al. 20081 



iMillion and AllenI [20091) . However, a recent search bv lAiello et aLr(|2009i ) of 
hard X-ray excess above the thermal emission in a large sample of clusters 
with Swift Burst Alert Telescope does not show significant nonthermal hard 
X-ray emiss i on. T he only exception is Perseus whose high-energy emission 
lAiello et afl ([2009D attributed to the central galaxy NGC 1275. In the next 



years, as more cluster outskirts are being observed with Suzaku, more con- 
clusive results on shocks and nonthermal components are expected. 

Compared to X-ray observations, the outer regions of clusters can poten- 
tially be better studied by Sunyaev-Zel'dovich (SZ) observations because the 
SZ effect depends only on the electron density to the first power, while the 
X-ray emission depends on the electron density squared. Facilities such as 
the Atacama Cosmology Telescope (ACT), the South Pole Telescope (SPT) 
and the Planck satellite are searching for the SZ signal of galaxy clusters and 
have found some interesting first results. For example, it was found that the 
SZ signal caused by clusters is by a factor of ~ 2 smaller than predicted by 
models of clu sters that suggest th at the pressure by the electrons has been 
overpredicted iLueker et al.1 ([2010[ ). This is presumably caused by a substan- 



tial nonthermal pressur e at cluster outsk irts, most of it is likely to be tur- 
bulent press ure Csee e.g.lShaw et all 2010l). This picture is also suppo rted by 



simulations ([Vazza et alJl2010llPaul et al.ll2010l iNaeai and Laull201lh which 



we review in Sec. 5. 

Given the advancements in the X-ray and SZ observations of the cluster 
outer regions, as well as the growing evidence of missing thermal energy in 
the ICM and the possible negative implications for cosmological tests, a more 
detailed study of the kinetic processes in cluster envelopes is necessary. 



3.2 7-ray Observations 



Recently, lAckermann et all ( 2010t ) reported on the search for GeV emis- 



sion from clusters of galaxies using the Large Area Telescope on the Fermi 



Gamma-ray Space Telescope. Only upper limits on the photon flux were re- 
ported in the range 0.2-100 GeV toward a sample of around 30 observed 
clusters (typical values (1 — 5) x 10~^ ph cm~^ s~^), considering both point- 
like and spatially resolved models for the high-energy emission. The authors 
concluded that the volume-averaged relativistic hadron-to-thermal energy 
density ratio is below 5% — 10% in several clusters. Also, us ing High Energy 
Stereo scopic System (H.E.S.S.) observations of Coma cluster. lAharonian et all 
()20Q9[ ) placed an upper limit ~ 10"^'^ ph cm~^ s~^ at photon energies above 
5 TeV for the Coma core of 0.2 degrees in radius constraining the multi- 
TeV particle population in the cluster. In summary, at this time the 7-ray 
observations remain inconclusive with respect to the origin of diffuse radio 
emission in clusters, but this is bound to change soon as the upper limits are 
expected to decrease. 



4 Cluster shocks in radio emission 

A number of diffuse, steep-spectrum radio sources without optical identifica- 
tion have been observed in galaxy clusters. The emission from these sources 
is synchrotron radiation, which indicates the presence of highly relativistic 
electrons and /iC range magnetic fields. Especially radio relics are thought to 
trace cosmological sho cks. Radio relics can be divided into two main groups 
(JKempner et al.ll2004l ). Radio gischt are large elongated, often Mpc-sized, ra- 



dio sources located in the periphery of merging clusters. They probably trace 
shock fronts in which particles are accelerated via the diffusive shock accelera- 
tion mechanism. Among them are do uble- relics with the two relics located on 
both s i des of a cluster center fe.g.. iBona fcdc ct al] 20091 van Weeren et al.l 



l2009bl IVeiituri et all [20071 iBagchi et al.ll200& iRottgering et al.lll997[ ). Ac 



cording to DSA, the integrated radio spectrum should follow a single power- 
law. Radio phoenices are related to radio-loud AGN. Fossil radio plasma from 
a previous episode of AGN activity is thought to be compressed by a merger 
shock wave which boosts, both, the magnetic field inside the plasma as well as 
the momenta of the relativistic particles. As a result the radio plasma bright- 
ens in synchrotron emission. In contrast to the radio gischt, the phoenices 
have a steep curved spectrurro indicating an aged population of electrons. 

The sizes of relics and the distance to the cluster centre vary signifi- 
cantly. Examples for radio relics with sizes of 1 Mpc or eve n larger have been 
observed in Coma (the prototype rel ic source 1253 -1- 275. (iGiovannini et al.l 
|1991[ ). Abell 2255 (iFeretti et al.ll997^ and Abell 2256 (IRottgering et al.lll994[ ). 



which contain both a relic and a halo (as do Abell 225, Abell 521, Abell 754, 
Abell 1300, Abell 2255 and Abell 2744). The cluster Abell 3667 contains two 
very luminous, alm ost symmetric relics with a separ ation of more than 5 Mpc 



(IRottgering et al.|[r9 97). as does ZwCl 2341.1-^0000 (Ivan Weeren et al.l2009J. 
IGiovannini et al.l20 101. and Abell 2345 and Abell 1240 (JBonafede et all2 00m. 



The clusters A115 and A1664 show re lics only at one side of the elongated 
X-ray distribution (JGovoni et al.l[200lr ). Recently, another double radio relic 
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Fi, (X v" , with a the spectral index 
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was found in the galaxy cluster ZwCl 0008.8+5215 (van Weeren et al. 2010) 
which we will discuss in detail below. 

The relic with the best evidence for shock acceleration found to date is lo- 
cated in the northern outskirts of the merging galaxy cluster CIZA J2242. 8+5301 
(z = 0.1921), see Fig. [51 The relic is located at a distance of 1.5 Mpc from the 
cluster center and spans 2 Mpc in length. The relic shows a clear unambigu- 
ous spectral index gradient towards the cluster center. The spectral index, 
measured between 2.3 and 0.61 GHz, steepens from —0.6 to —2.0 across the 
width of the narrow relic. The gradient is visible over the entire 2 Mpc length 
of the relic, something that has never been observed before. This is a crucial 
observation as it constitutes clear evidence for diffusive shock acceleration 
and spectral ageing of relativistic electrons in an outward moving shock. The 
relic is strongly polarized at the 50-60% level, indicating an ordered mag- 
netic field, and polarization magnetic field vectors are aligned with the relic 
(see Fig. [3]). In the southern part of the cluster a second fainter broader relic 
is found and the elongated radio relics are orientated perpendicular to the 
major axis of the cluster's elongated ICM. This is exactly as expected for 
a binary cluster merger event in which this second southern relic traces the 
shock wave that travels in the opposite direction from the first one. Fur- 
thermore, a faint halo of diffuse radio emission is seen extending all the way 
towards the cluster center connecting the two radio relics. This emission ex- 
tends over 3.1 Mpc, making it by far the largest known diffuse radio source 
in a cluster to date. 

The spectral index at the front of the relic is a = —0.6 ± 0.05. In simple 
shock acceleratio n theory, a is related to the Mach number via a = — (3 + 
M'^)/{2M'^ - 1) (JRosswog and Brugge"nll2007l ). which yields a Mach number 
of 4.6lQ'g. Subsequent hydrodynamic simulations of this shock front have 
indicated slightly lo wer Mach number of around 3. Using the Lx — T scaling 
relation for clusters ( Mar kevitchlll 9981 ) . we estimate the average temperature 



of the ICM to be '^ 9 keV. Behind the shock front the temperature is likely 
to be higher. Using the redshift, downstream velocity, spectral index, and 
characteristic synchrotron timescale, the width of the relic is given by 

^1/2^3/2 

;,enc«110kpcx-^— i^ (1) 

^ + -°CMB 

for a relic seen edge-on without any projection effects. i?cMB is the equivalent 
field strength of the IC scattering from the Cosmic Microwave Background, 
which is known, and thus the measurement of /relic from the radio maps, 
directly constrains the magnetic field. From the 610 MHz image (the image 
with the best signal to noise ratio and highest angular resolution), the relic 
has a deconvolved width (FWHM) of 55 kpc. Taking into account projection 
effects, van Weeren et al. (2010) conclude that the magnetic field strength at 
the location of the bright radio relic lies between 5 and 7 /iC As in A3667, 
this measurement points towards fiG - scale fields in merger shocks found at 
large distances from the cluster centre. 

Even though the sample of known radio relics is still small, one can start 
to find correlations between size, location and spectral index of these unique 
sources, which can eventually be compared to simulations of relic formation. 
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Fig. 2 Radio map (GMRT 325 MHz) of the relic in CIZA J2242. 8+5301. From 
Ivan Weeren et alJ ((MlOi) . 



The spectral index of the radio reUcs, versus the physical size is shown in 
Fig. 01 The projected distance from the cluster center is color coded. Van 
Weeren et al. (2010) find that on average the smaller relics have steeper 
spectra. Such a trend is in line with predictions from s hock statistics derived 
from cosmologica l simulations ( Skillman et al.l l2008l . iBattaglia et al.l l2009l 
iHoeft et al.ll2008l ). They find that larger shock waves occur mainly in lower- 
density regions and have larger Mach numbers, and consequently shallower 
spectra. Conversely, smaller shock waves are more likely to be found in cluster 
centers and have lower Mach numbers and steeper spectra. We note that 
more spectral index measurements of high quality are needed to confirm the 
correlation between physical size and spectral index. 

A further confirmation of the binary merger origin of most radio relics 
comes from measuring the angle of orientation of the relic's main axis with 
respect to the major axis of the X-ray surface brightness map of the host 
cluster. We have computed the angle a. between the major axis of the ICM 
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Fig. 3 Radio spectral index and polarization maps of the relic in CIZA 
J2242. 8-1-5301. Top: The spectral index was determined using matched observa- 
tions at 2.3, 1.7, 1.4, 1.2, and 0.61 GHz, fitting a power-law radio spectrum to the 
flux density measurements. Bottom: The polarization electric field vector map was 
obtained with the VLA at a frequency of 4.9 GHz. The length of the vectors is pro- 
portional the polarization fraction, which is the ratio between the total intensity 
and total polarized inten sity. A reference vector fo r 100% polarisation is drawn in 
the top left corner. From Fvan Weeren et al.l (|201CII 'l. 



and the line cluster centerradio relic center. The resulting histogram is shown 
in Fig. 5. From this histogram we see that relics are preferably found along 
the major axis of the ICM. This is in line with the simple picture that shock 
waves propagate outwards along the merger axis. We also calculate the angle 
(3 between the major axis of the ICM and the relics major axis. Fig. 18 shows 
that most relics are oriented perpendicular to the ICM major axis, also in 
agreement with a shock origin for radio relics. 

In the next years, simulations should be able to reproduce these trends, 
and thus constrain the efficiencies of shock acceleration and the magnetic 
field evolution in clusters. At the sa me time, it is pred icted that LOFAR will 
find of around 100 new radio relics ( Hoeft et al.ll2008l) . This will help address 
some remaining puzzles that surround radio relics, most importantly: 



Which processes accelerate electrons so efficiently at relatively low Mach 
number (M ~ 2 — 4) shocks? 
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Fig. 4 Spectral index of radio relics versus their size. Squares are the proposed 
radio phoenices. Diamonds represent the radio relics tracing merger shocks where 
particles are being accelerated by shock acceleration. The color coding is according 
to the projected distance from the cluster center. For the relics represented by black 
symbols w e could not obtain a reliable projected distance to the cluster center. From 
Ivan Wccren et all ((20098,). 




Fig. 5 Left: Schematic illustration of the angle between the major axis of the 
ICM and the line relic centercluster center (a), and the angle between the relic 
orientation and major axis of the ICM (/3). Top Right: Histogram of angles between 
the major axis of the X-ray emission and the line connecting the cluster center with 
the center of the relic. Bottom Right: Histogram of angle s between the major a xis 
of the X-ray emission and the relics major axis. From Ivan Weeren et al.l (| 20 111 ). 
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— What produces the magnetic fields inside relics? Both the inferred field 
strengths as well as the observed polarisation of the radio emission need 
to be explained. 

— Why do some relics have very sharp edges while others appear very fuzzy? 

— Under which conditions do relics form? When do we see single and when 
double relics? 

— Some relics appear to be connected to cluster-wide radio halo emission. 
Is there a physical connection between the two? 

In the next section, we will review the current state of cosmological simula- 
tions that include a treatment for shock acceleration of cosmic rays. 



5 What w^e know from simulations 

5.1 Shocks in cosmological simulations 

The statistics of cosmological shocks in the large-scale structure of the Uni- 
verse has been derived from s i mulations using, both , Eulerian hydrodynamic 

codes (e.g. lMiniat e t al."2000'.'200l'.'Rv u et al.l2003llKang et al.l2007liskillman erall 



[2008, Vaz za et al.l EoiO. 2009b.a. Molnar et al.ll2009D and smoo thed particle 



hydrodynamic codes (iPfrommer et a l. 2007, Hoeft et al. 2008) as we ll as a 



number of semi- analytical ( Gabici and Blasl.2003, .Keshet et al]|2003[ ) stud- 



ies. 

Ryu et al. (2003) performed uniform grid simulations with limited reso- 
lution, in which cluster core regions were not properly resolved, while sur- 
rounding outskirts were fairly well reproduced (see Fig. 6). Hence, the shocks 
identified in these simulations are mostly those in cluster outskirts. A cell is 
determined to have a shock if it meets the following requirements: 1. V-v < 0, 
2. VT -Vs > 0, 3. T2 > Ti and 4. p2 > pi- where v is the velocity field, T 
is the temperature, p is the density, and s = T/p"'^^ is, in X-ray astronomy 
parlance, the entropy. The subscripts 1 and 2 denote up- and downstream 
quantities, respectively. To find shocks, one loops through rows of cells along 
each of the coordinate axes and identifies one-dimensional shock structures in 
each direction. Fig. [7] shows the surface area, S, of identified shocks with the 
preshock gas temperature Ti > 10^ K, normalized by the entire volume of 
the simulation box, at the present epoch as a function of shock Mach number. 
The quantity S provides a measure of shock frequency. The Mach number 
of shocks expected to be found in cluster outskirts is low with M < 3. The 
freq uency incre a ses to w eakest possible shocks with M ~ 1. 

IVazza et al.l ( 2009H ) identifies shocks in cosmological simulations using 



a method based on velocities (instead of temperatures). They found that 
the overall differential distribution of shocks with their Mach number in 
the cosmic volume is very steep, with (3 ^ —1.6 (with MdN/dM ex M^ 
), and the bulk of the detected shocks at any Mach number is found in low- 
density regions. The Mach number distribution of detected shocks becomes 
increasingly steeper with ambient density: /? ~ —3 to —4 is found in clusters 
and their outskirts. They find that before the epoch of re-ionization, z > 
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Fig. 6 Two-dimensional slice around a region containing two clusters/groups at « = 
0. Top panels show the distribution of gas density (left) and temperature (right). 
Bottom panels show the locations of the shocks, which are color-coded according 
to shock Mach number M (left) or shock speed Vs (right). In the left panel, the 
colors are coded as follows: black for M > 100, blue for 30 < Af < 100, green for 
10 < M < 30, red for 3 < M < 10, and magenta for M < 3. In the right panel, the 
colors are coded as follows: black for Vs < 15 km s~^, blue for 15 < Vs < 65 km a~^, 
green for 65 < Vs < 250 km s~^, red for 250 < Vs < 1000 km s~^, and magenta for 
Vs < 1000 km s"\ 



6, roughly 30 % of the simulated volume is shocked. Then as soon as re- 
ionization plays a role, the temperature of the simulated volume increases 
and the Mach number distribution of shocks at redshift z ■^ 3 — 6 undergoes 
a dramatic change becoming very steep and dominated by weak shocks. The 
bulk of the energy in their simulations is dissipated in galaxy clusters which 
contribute about 75 % of the total energy dissipation (about 80 % if the 
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Fig. 7 Frequency of shocks expected to be found at z = in cluster outskirts as a 
function of M. Figure from D. Ryu. 



contribution from cluster outskirts is included), while filaments contribute 
about 15 % of the total energy dissipation. 

In line with previous numerical studies, relatively weak shocks are found 
to dominate the process of energy dissipation in the simulated cosmic volume, 
although they find a larger ratio between weak and strong shocks with respect 
to previous studies. The bulk of energy is dissipated at shocks with Mach 
number M r-^ 2 and the fraction of strong shocks decreases with increasing 
density of the cosmic environments. 

ISkillman et al.l (|2008l) improved the method employed by Ryu et al. (2003) 
that can produce errors when examining shocks whose direction of motion 
is not oriented along a coordinate axis. They find that the Mach number 
evolution can be interpreted as a method to visualize large-scale structure 
formation. Shocks with Af < 5 typically trace mergers and complex flows, 
while 5 < Af < 20 and Af > 20 generally follow accretion onto filaments and 
galaxy clusters, respectively. 
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h"'' Mpc 

Fig. 8 Radio emission from a galaxy cluster taken from the Mare Nostrum sim- 
ulation and resimulated at higher resolution. The vectors denote the degree and 
the direction of polarisation of the radio emission. The contours indicate the X-ray 
emission. (From M. Hoeft) 



5.2 Cosmological simulations with particle acceleration 



In recent years a number of cosmological simulations have started to include 
some treatment of cosmic ray physics with particle acceleration at shocks. 

Using a sho ck-finding formalism tuned to smoothed-particle hydrodynam- 
ics simulations. iHoeft et al.l ( 2008t ) have analyzed the MareNostrum Universe 
simulation which has 2 x 1024'^ particles in a 500 h~^ Mpc box. In a ddition, 
they have used the formalism derived in iHoeft and Briiggenl ( 20071 ) to pro- 
duce artificial radio maps of massive clusters. Several clusters were found to 
show radio objects with similar morphologies to observed large-scale radio 
relics (see Fig. [5]) , whereas about half of the clusters show only very little ra- 
dio emission. In agreement with observational findings, the maximum diffuse 
radio emission of galaxy clusters depends strongly on their X-ray tempera- 
ture. Moreover, it was found that the so-called accretion shocks cause only 
very little radio emission. 

Using a particle-mesh and Eulerian hydrodynamics code, Ryu et al. have 
been studying the gas thermalization and CR acceleration efficiencies of 
shocks in cosmological grid simulations. In order to quantify the energy 
dissipation at cosmological shocks, the incident shock kinetic energy flux. 
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-Fkin — (l/2)pivf where pi is the preshock gas density, and Vs the shock 
velocity, was calculated for each shock, and the results are shown in Fig. [T] 
Then, the kinetic energy flux through shock surfaces, normalized by the 
entire volume of the simulation box, as a function of shock Mach number 
was calculated. The resulting kinetic energy flux is shown by the dotted 
line in Fig. [HI The kinetic energy flux through shock surfaces is larger for 
weaker shocks; energetically weaker shocks are more important processing 
more shock energy. The acceleration efficiency for cosmic rays, 77, depends 
on the Mach number and the injection parameter, es, which is defined as 
eB = Bo/B±, the ratio of the mean magnetic field strength aligned with 
the shock normal to the amplitude of the postshock turbulent wave field. In 
addition, it is assumed that the CR population is isotropized with respect to 
the local Alfvenic wave turbulence, which would in general drift upstream at 
the Alfven speed with respect to the bulk plasma. This reduces the velocity 
difference between upstream and downstream scattering centers compared 
to the bulk flow, leading to less efficient shock acceleration. Moreover, the 
dissipation of Alfven turbulence heats the infiowing plasma in the precursor, 
which leads to weakening of the subshock strength. (See iKang and Jones! 
( 20071 ) for further discussions and references.) 



Fig. [To] shows the gas thermalization and CR acceleration efficiencies, de- 
fined as 6{M) = Fth/Fkin and r]{M) = -FcRZ-Pkin, respectively, where Fth 
and i^cR are the gas thermal and CR energy fiux, respectively, generated 
at shocks; they are from numerical simulations of quasi-parallel shocks with 
speed Vs = M ■ 150km s~^ (Ti = l O^K), and based on a shock accelera- 
tion model of iKang and Jonea (|2007f ) with e_B — 0.25. At strong shocks with 
M J> 10, the injection rate is high enough so that the CR acceleration ef- 
ficiency nearly saturates and becomes almost independent of the parameter 
eb- At weak shocks, on the other hand, the level of injection lies in the regime 
where the CR acceleration efficiency increases with the injection rate; conse- 
quently, 77 depends sensitivel y on cb ■ Not e that the value of rj presented here 
is ^ 1/2 of that presented in lKang et al.l (2007) for M <^ 5, while it remains 
about the same for stronger shocks. 

By adopting the efficiencies in Fig. [TUJ the gas thermal and CR energy 
fluxes dissipated at cosmological shocks as a function of shock Mach number 
were calculated in the same way the kinetic energy flux through shock sur- 
faces was calculated. The resulting thermal and CR energy fluxes are shown 
with dashed and solid lines in Fig. [HI For the generation of thermal energy, 
still weaker shocks are more important; shocks with M <^ 2.5 contribute the 
most to the thermal energy. On the other hand, the generation of CR en- 
ergy is more efficient at shocks with M Jj 2. These indicate that in cluster 
outskirts, while shocks with M < 2 are more abundant, shocks with M > 2 
would be more frequently detected as radio relics. We note that the results 
in Fig. ini are consistent with those estimated for cosmological shocks in all 
the regions of the IGM; while shocks around M ~ 2 contribute most to the 
generation of gas thermal en ergy, shocks around M ^ 3 contrib ute most to 
the generation of CR energy ()Rvu et al.ll2003l |KanieraL||2003). 

At weak shocks, DSA is known to be rather inefficient and the CR pressure 
remains dynamically insignicant, partly because the injection from thermal 



19 







t3 



UO 

o 



-4 



n I I I I I I I I I I r 



n I I r 



I I 




3.5 



Fig. 9 Dotted line: Kinetic energy flux per comoving volume passing through sur- 
faces of shocks in cluster outskirts as a function of M at z — 0. Dashed line: Gas 
thermal energy dissipated at the shocks as a function of M. Solid line: CR energy 
accelerated at the shocks as a function of M. All the quantities are in units of lO'*" 
ergs s"^(/i~^Mpc)"^ 



to nonthermal particles is inefficient (e.g., Kang et al. 2002). Recently, Kang 
& Ryu (2010) found that at weak; shocks expected to form in clusters much 
less than 10"'^ of particles are injected into CRs and much less than ~ 1% 
of the shock ram pressure is converted into the downstream pressure of CR 
protons, so the particle acceleration is virtually negligible. 

This has been followed up in recent AMR simulations by Vazza (2011 in 
prep.) that implement prescriptions for CR acceleration at shocks. This work 
suggests that the shock acceleration efficie ncy needs to be incr eased in order 
to explain the occurrence of radio relics. iKang et al.l ( 2007t ) have already 
suggested that that pre-existing populations of cosmic rays can boost the 
acceleration efficiency significantly, and this may be part of the answer. Thus 



20 



0.1 



0.01 



0.001 



— 


' 


1 


' 




1 ' 1 




— 


- 




/^.e-- 


--0-- 


-e- 


e-- 


--<)-- 


Z^ 




£f 








^•r^ 






/ 


^ 






^ 


^7] 




— 


- / 








Mi 






— 








/' 








- 


- 




/ 1 




1 


1 , 1 




- 



4 



6 



8 10 



M 



Fig. 10 Gas thermalization efficiency, 5(M), and CR acceleration efficiency, ri{M), 
as a function of Mach number. Open and filled circles are the values estimated from 
numerical simulations based on a shock acceleration model and dashed and solid 
lines are the fits. 



the observations of cluster radio relics can inform us about the microphysics 
in weak shocks. This will be the subject of the next section. 



6 Microphysics of shocks 

6.1 Plasma physics in CR- modified shocks 



The heating processes in collisionless shocks are complex phenomena. The 
irreversible transformation of a part of the kinetic energy of the ordered 
bulk motion of the upstream flow into the energy of the random motions 
of plasma particles in the downstream flow in collisionless shocks is due to 
the plasma instabilities in the thin dissipation region where the wave-particle 
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interactions provide momentum and energy redistribution between different 
components. 

The front of a strong coUisionless shock wave may consist of a precursor 
and a viscous velocity discontinuity (subshock) of a local Mach number that 
is smaller than the total Mach number of the shock wave (Fig. [TT|) . The 
compression of matter at the subshock can be much lower than the total 
compression of the medium in the shock wave allowing for a high compression 
in the precursor. We refer to such shocks as CR-modified. The standard 
collisional model of a strong single-fluid shock predicts a particle temperature 
kT -^ (3/16) -TOUgj^ for 7^ = 5/3. However, this is not generally valid in multi- 
fluid plasma and CR-modified shocks since the total shock compression ratio, 
Rt{vsh), depends on the shock velocity, Vgh, in the strong shock limit. 

In a simplified steady-state model, a strong CR-modified shock can be 

parameterized by the total Mach number of the shock A^tot and the Mach 

(2) 
number of the subshock A^sub- Then the downstream ion temperature T; 

can be estimated for the CR modified shock of a given velocity Ugh if the 

total compression ratio, Rt, is known: 

7g-Kt ("shj (7g - Ij-^sub + 2 

The total compression ratio Rt in Eq. [2] depends on the precursor heating 
rate and it is the main parameter in the steady state CR modified shock 
to det ermine the postshock ion temperature T; (|Bvkovll2005l iBvkov et al.l 

l2nn8aD . 

Shocks may also generate turbulence through cosmic-ray driven instabil- 
ities. In the shock precursor, this turbulence is certain to play a critical role 
in non-linear models of strong, CR-modified shocks. Different non-linear ap- 
proaches to the modeling of the larg e-scale structure of a shock undcrgoin. 



efficient cosmic ray acceleration (^e.g..lBell and Luceldl2 001. Am ato and Bias: 



2006llVladimirov et al.ll2006llBeUll2004IVladimirov et al.„2008K,Zirakashvih et al.l 



2008L iKang et al.l 120091 ) have been proposed. Many of these have predicted 
the presence of strong MHD turbulence in the shock precursor. An exact 
modeling of the shock structure in a turbulent medium, including nonthermal 
particle injection and acceleration, requires a nonperturbative, self-consistent 
description of a multi-component and multi-scale system including the strong 
MHD-turbulence dynamics. 

A Monte-Carlo model of nonlinear diffusive shock acceleration accounting 
for magnetic field amplification through resonant instabilities induced by 
accelerated particles, and including the effects of dissipation of turbulence 
up stream of a shock and the subsequent precursor plasma heating was studied 
by IVladimirov et al.l ()2008l ). Feedback effects between the plasma heating 
due to turbulence dissipation and particle injection are strong, adding to the 
nonlinear nature of efficient shock acceleration. 

De scribing the turbulence damping in a parametrized wav. lVladimirov et al.l 
( 20081 ) reach two important results: first, for conditions typical of super- 



nova remnant shocks, even a small amount of dissipated turbulence energy 
(~ 10%) is sufficient to significantly heat the precursor plasma, and second, 
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the heating upstream of the shock leads to an increase in the injection of ther- 
mal particles at the subshock by a factor of a few. In their results, the response 
of the non-linear shock structure to the boost in particle injection prevented 
the efficiency of particle acceleration and magnetic field amplification from 
increasing. The non-linear models models with turbulence generation and 
dissipation may lead to a scenario in which particle injection boost due to 
turbulence dissipation results in more efficient acceleration and even stronger 
amplified magnetic fields than without dissipation. In Fig. [TT] the effects are 
illustrated for a stationary, non-linear Monte-Carlo simulation of a forward 
shock in a supernova remnant of velocity of 5,000 km s~^ accelerating par - 
ticles up to maximal energies of about 10^ GeV ( Vladimirov et ah! |2008() . 



The shock velocity is higher than is expected for cluster accretion shocks 
but the simulation can be used to investigate the qualitative behavior and to 
obtain the appropriate scalings. The maximal energy depends on the ratio 
of the shock scale size and particle escape physics that is simply modeled 
by a free escape boundary in the Monte Carlo simulation. Ultra-high energy 
CR acceleration abo ve 10 ^ GeV by cluster accretion shocks was shown by 
iNorman et al.l (1993) and iKang et al.l ( 19961 ) to be a plausible scenario. In 



that case the scale size of CR precursor with magnetic field amplification can 
be about 100 kpc. 

IVladimirov et al] (|2008f) demonstrated that the effect of turbulent dissi- 
pation on the thermal plasma is evident in the values of the pre-subshock 
temperature Ti, the downstream temperature T2, and the volume-averaged 
precursor temperature (r{x < 0)) (the averaging takes place over the CR pre- 
cursor region). The temperatures were calculated from the thermal particle 
pressure using the ideal gas law. The value of the pre-subshock temperature 
Ti depends drastically on the level of the turbulent dissipation an, increas- 
ing from an = to an — 0.5 by a factor of 11 in the case of T = 10® K. It 
is less in the case of T = 10® K than for T ^ 10"* K because the efficiency of 
the CR streaming instability in generating the magnetic turbulence is less for 
the smaller Mach number. The values of the temperature as high as Ti are 
achieved upstream only near the subshock; the volume-averaged upstream 
temperature, {T{x < 0)), is significantly lower. The factor by which the av- 
erage temperature {T(x < 0)) increases in case of T ■^ 10® K is about 2.3. The 
multi-fluid processes described above can preheat the accreting gas. An im- 
portant prediction of the models of strong shocks is the possibility to amplify 
an initial seed magnetic field by a few orders of magnitude. CR currents and 
CR pressure gradients upstream of the strong shock can drive magnetic fluc- 
tuations on the shock precursor scale length. The precursor scale size is 10^ 
times larger than the subshock transition region where strong, small-scale 
magnetic field fluctuations are directly produced by instabilities of super- 
Alfvenic bulk plasma flows. These small-scale fluctuations are responsible for 
bulk dissipation and the adiabatic amplification of the transverse magnetic 
field in collisionless shocks. From Fig. [TTl it is clear that the magnetic field 
amplification factor depends on the non-linear processes of particle injection 
through turbulent dissipation in the shock precursor. 

Above, we discussed the ion temperature in multi-fluid shocks. However, 
cluster X-ray spectra depend primarily on the electron temperature, while 
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Fig. 11 Results of non-linear simulations in the case of the far upstream tempera- 
ture T = 10® K with different values of turbulence dissipation parameter an taken 
from Vladimirov et al. (2008). The solid, dashed and dotted lines correspond, re- 
spectively, to Oh ~ 0, 0.5 and 1.0. The plotted quantities are the bulk flow speed 
u{x), the effective amplified magnetic field Bcsix) and the thermal gas tempera- 
ture T{x). The shock is located at a; = 0. There is a change from the logarithmic 
to the linear scale at x = —0.05. The distances are measured in relativistic proton 
gyroradius unit r^o- 



MHD-type wave dissipation in the shock precursor may preferentially heat 
ions. Hovifever, the complete e — i Coulomb equilibration requires the sys- 
tem age Tci > 10^° Tq' /n, where the postshock density n is measured in 
cm~'^, the ion temperature Tq is in 10^ K, and Tci is in seconds iMewd (|l99Cl| ). 
Since the number density around the cluster virial radius is typically above 
10~^, cm~'^, the temperature equilibration scale length in the accretion shock 
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precursor is below 100 kpc. The electron temper ature downstream of the 
shock depends on collisionless electron hea ting e.g. (JBvkov and Uvarovlll999l 



iRakowski et al.l [20081 iBvkov et al.ll2008bl ). Non- resonant interactions of the 
electrons with strong nonlinear fluctuations generated by kinetic instabilities 
of the ions in the transition region inside the shock front may play the main 
role in the h eating and pre-acceleratio n of the electrons, as it was shown in 
the model bv lBvkov and Uvarovl (|l999l) . They calculated the electron energy 



spectrum in the vicinity of shock waves and showed that the heating and 
pre-acceleration of the electrons occur on a scale of the order of several hun- 
dred ion inertial lengths h in the vicinity of a viscous discontinuity. Although 
the electron distribution function is significantly out of equilibrium near the 
shock front, its low energy part can be approximated by a Maxwellian dis- 
tribution. 

The effect of the non-adiabatic electron heati ng efficiency, /3e, on the de - 
gree of non-equipartition was studied recently by IWong and SarazinI (|2009l) . 



They have shown that for a cluster with a mass of Mvir ^ 1.2 x IO^^Mq , 
electron and ion temperatures differ by less than a percent within the virial 
radius. The difference is 20% for a non-adiabatic electron heating efficiency 
of 1/1800 to 0.5 at ~ 1.4 of the virial radius. Beyon d this radius, the non- 
equipartition effect depends rather strongly on (3^ (Wong and Sarazin 2009) , 
and such a strong dependence at the shock radius can be used to distinguish 
shock heating models or constrain the shock heating efficiency of electrons. 

A s opposed to semi-analytic kinetic theory methods ()Achterberg et al.l 

I2OOI ^ Keshet and Waxman _ 2005^ or Monte Carlo test -particle simulations 



(JNiemiec and Ostrowskj20M . lEllison and Doublel2004[ ) . PIC simulations can 
tackle the problem of particle acceleration in shocks without the need for sim- 
plifying assumptions about the nature of the magnetic turbulence or the de- 
tails of wave-particle interactions. Multi-dimensional PIC simulations o f rela - 
tivistic un-magnetized shocks have been p r esent ed, e.g. bv lSpitkovskvl (120081 ) 



tivistic un-magnetized snocKS nave been p r esent ea, e.g. pv lbpitkovskvl ll/UUal l 
for electron-ion flows. iKato and Takabg ()2010[ ) present a two-dimensional 
particle-in-cell simulation to investigate weakly magnetized perpendicular 
shocks. The simulated thermal energies of electrons and ions in the down- 
stream region are not in equipartition but their temperature ratio /3e = 
Te/Ti ~ 0.3 - 0.4 is high enough indicating rather efficient non-adiabatic 
heating of the electrons in the shock transition region. 

Using 2.5D PIC simulations, ISironi and Spitkovskvl ([2010) find that in 
subluminal shocks, where relativistic particles can escape ahead of the shock 
along the magnetic field lines, ions are efficiently accelerated via a Fermi- 
like mechanism. The scattering is provided by short- wavelength non-resonant 
modes produced by Bell's instability (Bell 2004), whose growth is seeded by 
the current of shock-accelerated ions that propagate ahead of the shock. Up- 
stream electrons enter the shock with lower energy than ions, so they are 
more strongly tied to the field. As a result, only about 1% of the incom- 
ing electrons are Fermi-accelerated at the shock before they are advected 
downstream, where they populate a steep power-law tail. Thus, efficient elec- 
tron heating is the universal property of relativistic electron-ion shocks, but 
signicant nonthermal acceleration of electrons is hard to achieve in magne- 
tized flows and requires weakly magnetized shocks, where magnetic fields 
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Fig. 12 Normalized post-shock gas entropy, Kjns — KqT j p-'''- ^', as a function of 
the CR particle energy escape flux Qeac carried by energetic particles. The dimen- 
sionless flux is defined as r^osc = Qosc/-?7kin, where TTkin ~ piv^^/2. The upper curve 
(dotted) corresponds to an effective adiabatic exponent 7g = 4/3 (relativistic gas), 
while the lower (solid) curve corresponds to 7g = 5/3. The postshock entropy is 
normalized to Kmiirjcsc = 0). 



self-generated via the Weibel instability are stronger than the background 
field. 



6.2 Entropy in CR-modifled Accretion Shocks 



A distinctive feature of CR-modified shocks is their high gas compression 
factor Rt{vsh) that can be much higher than the single fluid shock limit of 
(7g -I- l)/(7g — 1). Since the postshock gas entropy for a strong multi-fluid 
shock scales as Rt (wsh)"'''*'''^'' , it can be significantly red uced compared to th e 
single-fluid adiabatic shock of the same velocity {e.g., iBvkov et al.ll2008al) . 
The effects are due to energetic particle acceleration and magnetic field gen- 
eration. Energetic particles can penetrate into the shock upstream region. 
They are coupled with the upstream gas through fluctuating magnetic fields 
(including the Alfven waves). Magnetic field dissipation provides gas pre- 
heating and entropy production in the shock precursor. Energetic particles 
deplete the momentum distribution in the gas resulting in a higher gas com- 
pression and a reduced temperature. In Fig. 12 the post-shock gas entropy 
is shown as a function of the fraction rjesc ■ The curves were calculated for a 
strong shock with no precursor heating i.e. an = 0. 

In galaxy clusters, entropy profiles have been measured by iPratt et al.l 
( 20101 ) at least out to i?iooo in 31 nearby galaxy clusters from the XMM- 
Newton cluster structure survey (REXCESS) and out to i?5oo in thirteen 
systems. The effects of CRs on the entropy production in shocks illustrated 
above can be of importance for accretion shocks that are expected to be 
located at larger radii. It will require deeper exposures of the next generation 
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telescopes to measure the entropy profiles at these radii. The entropy in the 
inner cluster regions can be also affected by internal shocks if the shocks are 
efficient accelerator s of CRs. The current entrop y profile simulations (e.g. 
iBorgani et al.l l2008l iBorgani and Kravtsovl l2009l) ought to be extended to 
include CR physics. 



7 Summary 

It is only now, with low-frequency radio telescopes, long exposures with high- 
resolution X-ray satellites and 7-ray telescopes, that we are beginning to 
learn about the physics of cluster outskirts. In the coming years SZ tele- 
scopes are going to deliver further insights into the plasma physics of these 
special regions in the Universe. The last years have already shown tremen- 
dous progress with detections of shocks, estimates of magnetic field strengths 
and constraints on the particle acceleration efficiency. The main points of this 
review are listed below: 

— Cosmological shocks are collisionless: Coulomb collisions are not sufficient 
to provide the viscous dissipation of the incoming ffow, and collective 
effects play a major role. 

— Collisionless shocks generate energetic, non-thermal particles that can 
penetrate far into the upstream flow. The particles decelerate the flow 
and preheat the gas. They can also efficiently generate strong fluctuating 
magnetic fields in the upstream region. This turbulence, which may result 
from cosmic-ray driven instabilities, in the shock precursor, is certain to 
play a critical role in non-linear models of strong, CR-modified shocks. 

— A few merger shocks have been identified in X-rays, exhibiting both a 
sharp gas density edge and an unambiguous temperature jum p: in the 
"bull et cluster" , IE 06 5756 (Markevitch et al. 2002) , A520 (Markevitch et al.l 
|2005^ and Abell 2146 (iRussell et al.ll2010D . 



On the basis of X-ray observations, microphysical properties of shocks 
can be derived. Markevitch et al. (2006) find that the temperatures across 
merger shocks are consistent with instant heating; equilibration between 
electrons and protons on the coUisional timescale is excluded. The electron 
temperature downstream of the shock depends on collisionless electron 
heatin g e.g. (JBvkov and Uvarovlll999L iRakowski et al.ll2008L iBvkov et al.l 



I2OO8H ) . Non-resonant interactions of the electrons with strong nonlinear 
fluctuations generated by kinetic instabilities of the ions in the transition 
region inside the shock front may play the main role in the heating and 
pre-acceleration of the electrons. Future X-ray observations extending to 
the virial radius or even close to the shock radius should be able to detect 
signatures of non-equipartition at shocks. 

Suzaku observations provide evidence for departures from hydrostatic 
equilibrium around and even b efore the viri al radius as it is seen in the 
galaxy cluster Abell 1795 {e.g., iBautz et al.|[2009^ . Evidence for nonther- 
mal pressure support suggests that bulk motions from mergers could be 
making a significant contribution to the gas energy in the outskirts of this 
cluster (|George et al.ii200&) . 
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Recently discovered radio relics provide very strong evidence for shock 
acceleration at merger shock waves in galaxy clusters. Van Weeren et 
al. (2010) find that on average smaller radio relics have steeper spec- 
tra. Such a trend is in line wit h predictions from shock statistics derived 
from co smological simulations ( Skillman et al.l l2008'.'Ba ttaglia et al.ll2009l 



iHoeft et al. 2008) . They find that larger shock waves occur mainly in 
lower-density regions and have larger Mach numbers, and consequently 
shallower spectra. 

— Observations of radio relics in combination with X-ray observations sug- 
gest magnetic field strengths in relics of > 3 /xG and efficiencies of 0.2% 
of the kinetic energy dissipated in a shock going into relativistic electrons 
(Finoguenov et al. 2009, van Weeren et al. 2010) . 

— Radio relics also suggest that the efficiency of shock acceleration at weak 
shocks is higher than predicted in theories of diffusive shock acceleration. 

— The bulk of the energy in their simulations is dissipated in galaxy clus- 
ters which contribute about 75 % of the total energy dissipation (about 
80 % the contribution from cluster outskirts is included) , while filaments 
contribute about 15 % of the total energy dissipation. In agreement with 
observational findings, the maximum diffuse radio emission of galaxy clus- 
ters depends strongly on their X-ray temperature. While shocks around 
M ~ 2 contribute most to the generation of gas thermal ener gy, shocks 
around M '-^ 3 contribu te most to the generation of CR energy (,Rvu et al.l 
l2003l[Kanget al.ll2007l) . 

— The effects of CRs on entropy production in shocks as illustrated above 
can be of importance for accretion shocks that are expected to be located 
at larger radii. It will require deeper exposures of the next generation 
telescopes to measure the entropy profiles at these radii. 

Many open questions remain. So far still little is known about shock acceler- 
ation and magnetic field generation in clusters, nor about the microstructure 
of shocks in very dilute plasmas. 

Most likely, all these processes are related and their understanding re- 
quires input from, both, observations and theory. The next years will see a 
rapid growth in observational data. On the theoretical side, we expect sub- 
stantial progress from new techniques that simulate particle acceleration in 
a magnetohydrodynamical setting. Fully kinetic particle-in-cell simulations 
provide a powerful tool for exploration of the structure of collisionless shocks 
from first principles, thus determining self-consi stently the interplay between 
shock-generated waves and accelerated particles (jSironi and Spitkovskvl2010l ). 



Results from these simulations can then serve as input for cosmological sim- 
ulations that include CR-physics. 
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